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> ■ ABSTRACT 

o 

1/'^ , Aims. To detect line effects using spectropolarimetry in order to find evidence of rotating disks and their respective 

■ symmetry axes in T Tauri stars. 
0^ ' Mettiods. We used the lAGPOL imaging polarimeter along with the Eucalyptus-IFU to obtain spectropolarimetric 

measurements of the T Tauri stars RY Tau (two epochs) and PX Vul (one epoch). Evidence of line effects showing a 
loop on the Q— U diagram favors a compact rather than an extended source for the line photons in a rotating disk. 
In addition, the polarization position angle (PA) obtained using the line effect can constrain the symmetry axis of the 
disk. 

Results. RY Tau shown a variable Ha double peak on 2004-2005 data. Polarization line effect is evident on the Q— U 
^ , diagram for both epochs confirming a clockwise rotating disk. A single loop is evident on 2004 changing to a linear 

■ excursion plus a loop on 2005. Interestingly, the intrinsic PA calculated using the line effect is consistent between our 
' two epochs (~167°). An alternative intrinsic PA computed from the interstellar polarization corrected continuum and 
I averaged between 2001-2005 yielded a PA~137°. This last value is closer to be perpendicular to the observed disk 

direction (^^25°) as expected by single scattering in an optically thin disk. For PX Vul, we detected spectral variability 
in Ha along with non- variable continuum polarization when compared with previous data. The Q— U diagram shows a 
well-defined loop in Ha associated to a counter-clockwise rotating disk. The symmetry axis inferred by the line effect 
^ ' has a PA'~91° (with an ambiguity of 90°). Our results confirm previous evidence that the emission line in T Tauri stars 

0^ , has its origin in a compact source scattered off a rotating accretion disk. 

O . 

, Key words, polarization - stars: individual; RY Tau, PX Vul - stars: circumstellar matter 

rn ^ 

1. Introduction nate from a compact source, the polarization pattern will 

OO depend on the specific geometry and bulk motions of the 

On the last years spectropolarimetry has become a powerful scattering particles surrounding the compact line source. In 

^ and unique tool for studymg envelopes (or disks) around ^ rotating disk-like configuration, the subsequent breaking 

• unresolved young stellar objects (YSOs). In this sense, the left-right reflection symmetry in the velocity fields leads 

^ analysis of the Stokes parameters Q- U diagram has shown ^ changing in polarization angle with wavelength, which 

. to be a fundamental diagnost ic too l for polari zation Ime appears as a loop on the (?-[/ diagram (Vink et al. 12001) . 

; detections (Oudma ijer fc D rewll9991 Vmk et al. 120021 12U03| interestingly, the sense of rotation of the envelope/disk as 

— hereinafter V03, |2005al). gggj^ from the Earth also can be retrieved of the direction in 

Considering the line effects, the most simple case hap- which the PA rotates arou nd the loop on the Q- [/diagram 

pens when the emission lines are formed over a much larger (Poeckert & Marlborough IT977|) 
volume than continuum. In this case, the circumstellar scat- 
tering material will polarize the continuum light more than 

the line photons. The last ones will add mainly unpolarized The survey of Vink et al. (i2005al) found that the most 

flux because they will be less scattered, and therefore the (9/^0) ^he T Tauri stars (TTS) shown hue effect associ- 

net polarization across the emission Hue will be reduced, ^^d with a compact source of line photons that is scattered 

This feature is usually name d the depolarization Hue ef- ^ rotating accretion disk. This means that loops on the 

feet (Clarke & McLean [19^ and a linear excursion is ex- Q~ U diagram must be a well-defined signature for TTS. 

pected on the Q- U diagram with the position angle (PA) The statistics clearly has to be increased to confirm this 

unchanged. issue. 

Models to explain the line effect including rotation also 

have been e xplored (Po eckert fc Marlborough [1977[ Wood In this letter, we report spectropolarimetric observa- 

et al. 119931 Vink et al. I2005bp . If the line photons origi- tions around Ha of the TTS RY Tau and PX Vul. These are 

the first results of an ongoing spectropolarimetric survey of 

Send offprint requests to: A. Pereyra, e-mail: pereyra@on.br YSOs. The observations and data reduction are presented 

* Based on observations obtained at the Observatorio do Pico in § [H The results are shown in § [S] A summary with the 

dos Dias, LNA/MCT, Itajuba, Brazil. final conclusions are drawn in § 01 
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2. Observations 

The observations were niade in two runs on 2004 and 2005 
using the 1.6m telescope at the Observatorio do Pico dos 
Bias (OPD), Brazil. We used lAGPOL, the lAG i maging 
polarimeter (Magalhaes et al. 119961 Pereyra I2000p . along 
with the Eucalyptus-IFU (EIFU, de Oliveira et al. [2003]) . 
EIFU is an integral field unit composed of an array of 
32 X 16 50 /im fibers that covers a field of 30" x 15" on sky 
with a scale of 0.93" per pixel (Fig. [Tl left). The detector 
used was a Marconi 2048 x 4608 pixels back-illuminated 
CCB with 13.5 /im^ per pixel. The used spectral range for 
EIFU was -600 A centered in Ha (i?=4000) that yields a 
resolution of —0.3 A per pixel. 

We used lAGPOL in linear polarization mode with 
a Savart plate as analyzer and an achromatic A/2- 
waveplate as a retarder. Each measurement consisted in 
four or eight waveplate positions separated by 22? 5. The 
lAGPOL+ElFU setup permits that the beam collected by 
the telescope and divided in two orthogonal polarization 
beams by lAGPOL will be projected on the EIFU fibers 
array. Then, spectropolarimetry can be done using the or- 
dinary and extraordinary beams that yield the -o and -e 
spectra in each waveplate position, respectively. A log of 
observations is shown in Table [T] 

We used standard ZR^jQ procedures for IFU reduc- 
tions in each waveplate position image including bias and 
flatfielding corrections along with wavelength calibrations. 
Then, a special routine was developed to extract and stack 
fibers for the -o and -e beams and the sky region (Fig. (TJ 
right). The optimum aperture radius for the fibers extrac- 
tion was selected by minimizing the polarimetric errors. 
Therefore, three spectra were constructed in each wave- 
plate position: 0(A), E(A) and sky(A) where special care 
was taken to remove cosmic rays. The sky spectrum was 
constructed stacking fibers away from the -o and -e beams 
(Fig. [1] right) and subtracted of the -o and -e spectra in 
each waveplate position image. Several relative positions 
for the sky were tested with similar results. After that, we 
used SPECPOLQ package to construct the flux, polariza- 
tion, polarization PA and polarized flux spectra for a proper 
binning. This package allows binning the spectra using a 
variable bin size set by a selected and constant polarization 
error per bin. The errors are obtained from the residuals 
of the observations at each waveplate position image with 
respect to the expected cosinusoid curve. In general, they 
are consistent with the photon statistics. 

Observations of polarized standard stars in each run 
yielded the position angle correction to the equatorial sys- 
tem (PAcorr) and permitted to check the polarization effi- 
ciency of our setup. Our measurements are in agreement 
with the literature values within one sigma. In addition, 
the instrumental polarization was checked using unpolar- 
ized stars and estimated to be lower than 0.2%, and there- 
fore no correction was needed. A summary of the calibra- 
tion data is shown in Table [H Typical signal-to- noise ratios 
(per bin) gathered by our setup for the sources in this work 
are between 170 and 250. We checked for the polarization 
bias effect and it shown to be negligible. 

^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 

^ the original version of SPECPOL was written by A. Carciofi. 



Table 1. Log of observations. 



Object 




Date 




Total IT'^ 




(mag) 






(s) 


T Tauri 


RY Tau 


10.2 


2004/sep/05 


8 


2400 






2005/sep/03 


4 


2400 


PX Vul 


11.7 


2005/sep/03 


8 


4800 


standards 


HD 147084'' 


4.6 


2004/sep/05 


8 


1600 


HD 187929'^ 


3.9 


2005/sep/03 


8 


640 


HD 10476'= 


5.2 


2005/sep/03 


8 


800 



(a) from Simbad; (b) number of waveplate positions; (c) total 
integration time (N x individual IT by waveplate position) ; (d) 
polarized; (e) unpolarized. 



Table 2. Calibration summary. 



Object Pob3° PAob,'^ PAiit" PAeorr 
(%) (°) (%) (°) (°) 

HD 147084 4.5 (0.2) O 4.50 (0.01) 3iO 4-33.0 

HD 187929 1.5 (0.2) 37.9 1.62 (0.05) 91.0 +53.1 

HD 10476 0.2 (0.2) 0.02 (0.01) 

errors in parenthesis, (a) mean values for a speetrum binned witli a polar- 
ization error (per bin) of 0.2%, excepts HD10476 (unbinned). The quoted 
polarization error is the averaged error per bin. (b) References: HD 147084, 
Tapia ff988;l, R filter; HD 187929, Bailey & Hough 119821 . R filter; HD 
10476, Bastien et al. (1988), V filter. 

EB ^T|p EE -^r 

' Extracted fibers (aperture radius: 2pix = 2") 

Fig. 1. Left: 2-B reconstructed image from EIFU-f lAGPOL 
showing the ordinary (OB) and extraordinary (EB) beams 
for point sources in one of the waveplate positions. The 
field size is 30"xl5" with an scale of 0.93" per pixel. Right: 
Extracted fibers (in magenta) used to construct the ord., 
ext. and sky spectra in each waveplate position for an ap- 
propriate aperture radius. In this example, r=2pix or 13 
stacked fibers. The central fiber is shown in green for each 
spectra. 

3. Results 

3.1. RY Tau 

RY Tau is a classical TTS (F8III type, Mora et al. IM?T|) 
with a high optical variability and possibly associated with 
variable circumstellar obscuration as in UX Ori type (Beck 
& Simon I200ip . Spectral variability also has been observed 
(Vrba et al. 119931 Petrov et al. 1999) without a clear corre- 
lation between the Ha line strength and the optical mag- 
nitudes. 

RY Tau shows a historical strong linear polarization 
variability (Hough et al. 1981, Bastien 119821 Bergner et 
al. 119871 V03 and references therein) . The proposed polar- 
izing mechanisms include rotational modulations by cool 
and hot spots on the stellar surface (Stassun & Wood 
I1999P and variable extinction by a dusty disk (Menard 
& Bastien |I992j). Recently, robust evidence of a bipo- 
lar jet (PAjet =115°) was found in RY Tau (St-Onge & 
Bastien I2008P using Ha— continuum imaging that can be 
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Table 3. Results. 



Object 


Date 


p a 

obs 




PAint 






(%) 


n 


n 


RY Tau 


2004/sep/05 


1.8 (0.3) 


11 (5) 


171 (11) 




2005/sep/03 


2.6 (0.2) 


16 (2) 


161'' (4) 


PX Vul 


2005/sep/03 


3.2 (0.3) 


33 (2) 


91" (25) 



errors in parenthesis, (a) mean values for a spectrum binned with 
a polarization error (per bin) of 0.3%. The quoted polarization 
error is the averaged error per bin. (b) rebinned at 0.15% per 
bin. (c) or 91+90. 



used to infer the direction of the unresolved disk and con- 
strain the polarimetric information. 

Table [3] shows the observed continuum polarization for 
RY Tau integrated along the full range of our spectra. 
A slight enhancement of the continuum polarization level 
is observed with the PA practically unchanged (~14°) in 
one year. Our values are comparable with the data from 
Vink et al. (|2005aP between 2001 and 2003. The spectropo- 
larimetry around Ha is indicated in Fig. [51 The spectral 
range of the top spectra is only of 45A for a better visu- 
alization of the line feature. The spectra are binned using 
a variable bin size with a constant polarization error (per 
bin) of 0.3% and 0.15% for the 2004 and 200^ data, re- 
spectively. These errors are the best compromise between 
a minimum error (per bin) and line effect detection. The 
lower plots are the Q—U diagrams showing the line effect. 

On 2004, the total flux of Ha appeared as a double- 
peaked profile with the central dip at the continuum level 
and the redder peak more intense. On the other hand, the 
2005 data again shown a similar double-peak but with the 
central dip just below the blue peak. The AVVSO database 
(Henden i2007p shown that RY Tau in average faded from V 
-10.0 mag on 2004/sep to -10.8 mag on 2005/sep. The rel- 
ative intensities of the Ha emission components remained 
practically constant between the two epochs and can be 
classified as II— B type (Reipurth, Pedrosa & Lago I199"6|) . 
We did not find evidence of an increment of Ha when the 
star is dimmed (Petrov et al. 119991 ). 

We detected the Ha polarization line effect for RY Tau 
using the Q—U diagrams (Fig. [21 bottom) around the cen- 
tral dip (between 6560-6565 A). This wavelength range is 
also indicated on each spectrum (Fig. [2l top) for a better 
comparison. On 2004, the line effect produces a loop al- 
ready reported by V03. Nevertheless, the 2005 data show a 
stronger and more complex pattern that resembles a mix- 
ing between a loop and a linear excursion. The integrated 
continuum polarization immediately before and after of the 
line effect is the same in each epoch. Interestingly, the di- 
rection in which the wavelength increases around the loop 
(indicated by the arrow) has the same sense in our two 
epochs (and also in V03). Therefore, this confirms that the 
disk of RY Tau rotates in a clockwise direction as seen by 
an observer at the Earth. 

In order to determine the intrinsic PA (PAjnt) we used 
a linear fit on the points that display the line effect on the 

^ The Stokes parameter solution with few waveplate positions 
(as in 2005 data) trends to underestimate the polarization er- 
ror (per bin). In order to detect the line effect, the 2005 data 
was rebinned at a level that guaranteed the best compromise 
between the line effect detection and the minimum dispersion of 
the continuum polarization. 



Q—U diagram. Formally, the sloop of the fit yields two 
possible solutions (PAjnt and PAjnt +90°). In general, this 
ambiguity is solved by the solution more consistent with the 
relative positions of the interstellar polarization (ISP) and 
the observed polarization on the Q— U diagram (Schulte- 
Ladbeck et al. [T994|) . 

Fig. [31 shows the observed continuum polarization and 
the PAjnt directions obtained from the line effect for our 
two epochs (Table [3]). The ISP towa rd RY Tau (2.8% at 
PA=26°) computed by Petrov et al. (|1999p is also shown. 
Analysing this diagram it seems clear that the solutions 
with U < for PAjnt are more likely. These solutions 
are indicated on the last column of Table [3l As we can 
see, the PAjnt values are practically the same (—167°) be- 
tween our two epochs considering the individual errors. As 
mentioned above, the same procedure applied by V03 on 
2001 data yielded a PAjnt -^146° (also shown in Fig. ^ 
that is consistent with our results within —20°. Therefore, 
the three existing measurements using the line effect in Ha 
spectropolarimetry yielded a mean PAjnt —160° with a vari- 
ation range of —15°. 

Assuming the ISP from Petrov et al. (|1999p and dis- 
counting this value of the continuum polarization (Table[3]), 
the intrinsic polarization is 1.6% at PA=134° and 1.0% at 
PA =149° o n our 2004 and 2005 data, respectively. Vink et 
al. (|2005ap reported continuum polarization around Ha in 
four epochs between 2001 and 2003. Therefore, consider- 
ing the ISP used above, the mean intrinsic polarization for 
RY Tau between 2001 and 2005 is 1.3% at PA=137°. The 
variation range of the PAjnt is —25° in this interval. 

Comparing the mean values for the PAjnt by the above 
two methods, we obtained that they are reasonably con- 
sistent. In particular, the ISP corrected continuum method 
yields a PAjnt closer to be perpendicular to the apparent 
disk direction (PAdisk=25°) derived by the jet (St-Onge & 
Bastien [2008p . This fact reinforces the idea of an optically 
thin disk in RY Tau if single scattering is the prevailing 
mechanism. It is interesting to note the apparent align- 
ment between the disk direction and the ISP (see Fig. [3]) 
suggests a weak interplay between the local magnetic field 
and its presumed collimator effect on the jet (De Colle & 

3.2. PX Vul 

PX Vul (LHa 483-41) was initially cataloged as an early 
type star with circumstellar shell and Ha in emission at 
a distance of 420 pc and belonging to the R association 
Vul R2 (Herbst et al. I1982j) . Its position is coincident with 
the edge of the filamentary cloud L778 near its northern 
NH3 core (Myers et al.^ESI- Mora et al. (|M?T|) classified 
PX V ul as an early TTS (F3 Ve) and Hernandez et al. 
I^UU^ as a Herbig Ae/Be star (F3±1.5). Never theless , the 
last inferred masses (<2M0) by Eisner et al. (|2005p and 
Manoj et al. (|2006p seem more consistent with a TTS. 

Optical photometric variability (A V —0.3 mag) was re- 
ported by Herbst et al. (|1982p . Nevertheless, Eiroa et al. 
pOOSp classified PX Vul as non- variable {AV <0.05 mag). 
This object shows a moderate veiling level (r/{ =0.8±0.4), 
an important rotational velocity wsini =78±11 km s^^, and 
a higher mass-accretion rate of 1.3x10^^ M0 yr^^ (Eisner 

et ai. mm . 
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RY Tau 2004set05 binl 0.3% var delta+33.0 



4 




RY Tau 2005sep03 [1-4] 0.15% var delta +53,1 



6540 



6545 6550 6555 6560 6565 6670 6675 6580 6586 

Mi) 



0.00 0.01 0.02 0.03 0.04 
Q{x100%) 



^ 




6540 6545 6550 6555 6560 6565 6570 6575 6580 6585 
MA) 




0.01 0.02 0.03 
Q(x100%) 



Fig. 2. Spectropolarimetry of RY Tau around Ha on 2004/sep/05 (left) and 2005/sep/03 (right). Top: The vertical lines 
indicate the region (6560-6565 A) where the line effect is evident. The spectra show the total flux (top), polarization 
(upper middle), polarization PA (lower middle), and polarization flux (bottom) binned using a variable bin size with a 
constant polarization error (per bin) of 0.3% and 0.15% for the 2004 and 2005 data, respectively. Bottom: Q— U diagram 
showing the line effect (black dots). The direction in which the wavelength increases is also indicated (arrow) along with 
the mean values for the blue (white square) and red (black square) continuum. The errors plotted around the Hne 
effect are assuming photon statistics directly. 



The optical polarization of PX Vul was not found to 
be variable by Oudmaijer et al. (|200ip . Their mean R 
broadband polarization (3.8% at PA=27°) is comparable 
to our observed polarization integrated along the full spec- 
trum (Table [3]) , confirming therefore the polarimetric non- 
variability. 

The spectropolarimetry of PX Vul on 2005 /sep is shown 
in Fig.lH The total flux of Ha appeared as a double-peaked 
profile with the central dip at the continuum level. In rel- 
ative fiuxes, the red peak is ~2.1 times more intense than 
the blue one. For comparison. Ha shows a triple-peaked 
profile on 2004/jun (Eisner et al. '2005) being the redder 
peak more intense than the other two (1.8 and 3.9 times, 
respectively). The bluer peak is not present in our data 
(1.3 year after) and spectral variability seems to be present 
on the Ha emission changing from 11— Bm to 111— B type 
(Reipurth, Pedrosa & Lago 1996). This variability is also 
supported by the ratio between the redder peak and the 
continuum that fails from 2.9 (on 2004) to 2.3 (on 2005) 
(see Fig. El). 

As in § 13.11 we detected the Ha polarization line ef- 
fect for PX Vul using the Q—U diagram. In this case, the 
main line effect is on the red peak (around three points 



on diagram between 6563-6568 A) and appears as a loop 
consistent with the detection of a rotating disk. For a bet- 
ter comparison, we also plotted the points around the blue 
peak (6555-6563 A) where another and subtle effect can be 
present. However, the strength of this effect is of the order 
of our errors and therefore this was not considered. 

The linear fit on the points that display the loop on 
the Q—U diagram yields two possible solutions for the 
PAint (91 and 91+90°, Tablc[3]). As previous information is 
not available, the ambiguity on PAi„t cannot be solved for 
PX Vul. Following the sense of the PA rotation on the Q~ U 
diagram, the disk in PX Vul rotates in a counter-clockwise 
direction as seen from the Earth. 

4. Summary and Conclusions 

We shown here the first results of an ongoing spectropolari- 
metric survey of YSOs gathered with EIFU-I-IAGPOL at 
LNA. All the line detections shown in this work are signif- 
icant at least on one of the Stokes parameters with a 95% 
confidence level. 

The classical TTS RY Tau was observed in 2004 and 
2005. Interestingly, the double-peaked Ha emission profile 
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0.00 0.01 0.02 0.03 0.04 
Q (x100 %) 

Fig. 3. Q—U diagram for RY Tau showing the intrinsic 
polarization PA (sohd hnes) from the Ha hne effect on 
three epochs (2001/decfrom Vink et al. l2003t and, 2004/sep 
and 2005/sep from this work). The continuum polarization 
around Ha also is shown in each case (black dots). The in- 
terstellar polarization (Petrov et al. iI99 9) is also indicated 
(white dot) along with the PA disk (dashed line, St-Onge 
& Bastien I2008P . The arrows indicate the reference system 
(Q',U') corrected by the ISP. 



(II-B type) changed in ~one year from a central dip at the 
continuum level to a dip just below the blue peak. RY Tau 
faded about ~0.8 mag between our observations with the 
blue/red peak ratio practically unchanged and without ev- 
idence of an anti-correlation between Ha and the photom- 
etry as suggested by previous works. 

Analysing the polarization spectra on the Q—U dia- 
gram, RY Tau shown a well-defined line effect passing from 
a loop on 2004 to a more complex pattern (loop-|-linear 
excursion) on 2005. Our observations are consistent with 
previous spectropolarimetric data from V03 where loop de- 
tection was reported. This feature is interpreted as a signa- 
ture of a compact source for Ha— line photons. The intrinsic 
polarization PA obtained directly from the line effect did 
not seem change between the interval of our observation 
(PAint ~167°). We confirm that the disk in RY Tau rotates 
clockwise. Considering the continuum polarization, a slight 
enhancement was observed on 2005 with the PA practi- 
cally unchanged (~14°). The PAjnt computed by the ISP 
corrected continuum method is more consistent with the 
presence of an optically thin disk in RY Tau. 

For PX Vul, the comparison with previous spectra 
shown that spectral variability is present on the Ha emis- 
sion passing from a II-Bm to II-B type in ~1.3 yr. 
Nevertheless, the continuum polarization seems to be non- 
variable. We reported the detection of a rotating disk from 
the loop around Ha line on the Q—U diagram. The intrinsic 
polarization has a PA~91° (with an ambiguity of 90°). An 
independent determination of the foreground polarization 
is necessary to better constrain the PAint- The sense of the 



PX Vul 2005sep03 0.3% var delta+53.1 




6540 6545 6550 6555 6560 6565 6570 6575 6580 6585 



A (A) 




0.00 0.01 0.02 

Q (xlOO %) 



Fig. 4. Spectropolarimetry of PV Vul around Ha on 
2005/sep/03. Top: The huge vertical lines indicate the full 
region (6555-6568 A) where the line effect is evident. The 
thin vertical line indicates the limit (6563 A) between the 
blue and red line effect. The spectra are as in Fig.[2]binned 
using a variable bin size with a constant polarization error 
(per bin) of 0.3%. Bottom: Q— U diagram showing the line 
effect (6555-6563 A in white dots, 6563-6568 A in black 
dots). The direction in which the wavelength increases is 
also indicated (arrow) along with the mean values for the 
blue (white square) and red (black square) continuum. The 
errors plotted around the line effect are assuming 
photon statistics directly. 

rotation of the PA on the loop shows that the direction of 
the disk rotation of PX Vul seems to be counter-clockwise. 

The finding of the signature of a loop in RY Tau and 
PX Vul confirms the trend of a compact source of line pho- 
tons that is scattered off a rotating accretion disk in TTS. 
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to tlic PCI/MCT/ON program). A.M.M. acknowledges support 
from FAPESP and CNPq. Polarimetry at lAG-USP is supported by 
FAPESP grant 01/12589-1. 
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